Excitation-assisted inelastic processes in trapped Bose-Einstein condensates 
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We find that inelastic coUisional processes in Bose- 
Einstein condensates induce local variations of the mean- 
field interparticle interaction and are accompanied by the cre- 
ation/annihilation of elementary excitations. The physical 
picture is demonstrated for the case of 3-body recombination 
in a trapped condensate. For a high trap barrier the produc- 
tion of high-energy trapped single-particle excitations results 
in a strong increase of the loss rate of atoms from the con- 
densate. 
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Since the discovery of Bose-Einstein condensation in 
trapped ultra-cold atomic gases [Q, inelastic coUisional 
processes in these systems attract a lot of attention. 
Three-body recombination and two-body spin relaxation 
limit achievable densities of trapped clouds and 
spin relaxation in atomic Cs even places severe limita- 
tions on achieving the regime of quantum degeneracy . 
Theoretical studies of 3-body recombination and 2-body 
spin relaxation in ultra-cold gases provide a valuable in- 
formation on the mechanisms and rates of these processes 
(see Q for review, and [^,0 for the earlier work in hy- 
drogen). These studies rely on a traditional approach of 
coUisional physics in gases: The decay rates are found on 
the basis of the calculated probability of inelastic transi- 
tion for a two- or three-body collision in vacuum. 

In this Letter we find that inelastic coUisional pro- 
cesses in Bose-condensed gases induce local variations of 
the mean-field interparticle interaction and can become 
"excitation-assisted". We demonstrate this phenomenon 
for the case of 3-body recombination. The particles pro- 
duced in the course of recombination have a high kinetic 
energy and immediately fly away from the point of re- 
combination. Therefore, each recombination event leads 
to an instantaneous local change of the mean field, i.e. to 
a change of the field acting on the surrounding particles. 
In this respect the recombination leads to "shaking" of 
the system: The time scale on which the Hamiltonian 
changes is so short that the wavefunction of the system 
remains unchanged and, hence, corresponds to a super- 
position of many eigenstates of a new Hamiltonian. For 
this reason the recombination event can be accompanied 
by the creation or annihilation of elementary excitations. 

From a general point of view, the concept of "shaking" 
(sudden perturbations) was formulated by Migdal and 
used for calculating the ionization of an atom in /3-decay 
(see, e.g., [0). Sudden perturbations in condensed sys- 



tems are accompanied by many-body collective effects. 
Dilute Bose-Einstein condensates are unique examples of 
gases, as the mean-field interparticle interaction makes 
their behavior in many aspects similar to that of a solid. 
To some extent, the picture of inelastic processes in these 
gases resembles, for example, the absorption of light by 
impurity particles in solids, where the transfer of impu- 
rities to an excited state leads to a sudden local change 
of the polarization of the medium and, hence, to the cre- 
ation/annihilation of phonons (see, e.g., [p2[). 

A distinct feature of "shaking" in a gaseous conden- 
sate concerns the back action of the shaking-created ex- 
citations on the condensate. This effect can be strongly 
pronounced in a trapped condensate. For a high trap 
barrier created high-energy single-particle excitations are 
still trapped and collide with condensate particles, re- 
moving them from the condensate. Every collision thus 
produces already two energetic (trapped) atoms which 
again collide with the condensate etc. As a result, one 
has a cascade production of non-condensed atoms out 
of the condensate. Despite a small probability of creat- 
ing high-energy excitations in the recombination process, 
this mechanism can significantly increase the total loss 
rate of Bose-condensed atoms. 

We will consider the most important channel of 3-body 
recombination in a trapped Bose-condensed gas, that is 
the recombination involving 3 condensate particles. For a 
single recombination event, the local change of the inter- 
particle interaction, SHi, can be obtained by considering 
the center of mass of the recombination-produced fast 
atom and molecule as a force center. Just before the 
event there are 3 atoms at the point , and each particle 
j of the sample interacts with this point via the potential 
3C/(rj — Ti). After the event this interaction is equal to 
zero, since the fast atom and molecule fiy away from the 
area of recombination. We will use a point approxima- 
tion for the interaction potential: U{r) = US{r), where 
U = 4:7rh^a/m, with a being the s-wave scattering length 
and m the atom mass. Then, assuming that the conden- 
sate density uq in the BEG spatial region greatly exceeds 
the density of non-condensed atoms, we obtain 

6H^^-3jd^r,i'\r,)Uir,-r,)^ir,)=-3Ui'\r,)^ir,).il) 

1/2 

The field operator of atoms, ^=^'o+5'', where 4*0 — tiq is 
the condensate wavefunction, and 4*' is a non-condensed 
part of the operator. This part refers to quasiparticle ex- 
citations, and we linearize SHi with respect to 5*'. Omit- 
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ting the term —SnoU, which is decoupled from the quasi- 
particle excitations, we obtain 



= -3C/*o(*'^(rO + *'(r.)). 



(2) 



The force center represents a sort of a "hole" in the 
coordinate space. It has the mass 3m and can undergo 
a translational motion. Therefore, the recombination- 
induced (sudden) change of the Hamiltonian Hq of the 
excitations can be written as 



SH= / drAHr^){SH^-{fi^/6m)A}^r,), 



(3) 



where $(r) is the field operator of the force center. The 
first term in the rhs of Eq. (||) is related to the interparti- 
cle interaction, and the second term to the motion of the 
force center. 

We assume that the kinetic energy of fast particles pro- 
duced in the recombination process greatly exceeds any 
other energy scale in the problem, and hence the creation 
or annihilation of excitations does not influence the en- 
ergy conservation law for the recombination. Then, ac- 
cording to the general theory of sudden perturbations 
PI , in each recombination event the probability of tran- 
sition of the excitation subsystem to a new state /, char- 
acterized by a different set of quantum numbers for the 
excitations, is given by Wij — The symbol {i\f) 

stands for the overlap integral between the wavefunction 
of the initial state i, which is an eigenstate of the Hamil- 
tonian Ho, and the wavefunction of the state / which 
is an eigenstate of the new Hamiltonian Hq + 5H . As 
Wif = 1, the creation/annihilation of excitations does 
not change the total recombination rate. 

In Thomas-Fermi condensates the most important is 
the creation of excitations with energies of order the 
chemical potential /i or larger (see below). These excita- 
tions are essentially quasiclassical, and their de Broglie 
wavelength is much smaller than the spatial size of the 
condensate. Hence, the probability of recombination ac- 
companied by the creation/annihilation of the excitations 
can be found in the local density approximation. In other 
words, as well as the recombination without production 
of excitations, this process occurs locally at a given point 
r characterized by local values of the chemical poten- 
tial fi and condensate density uq. Hence, one can use 
the Bogolyubov transformation for the spatially homo- 
geneous case and represent ^^j^* in terms of the cre- 
ation/annihilation operators Sj^, Sk of excitations charac- 
terized by momentum k; 

1 /2 

«''t(r,) + *'(rO = -^X^ (^^) (&I + 6-k)cxp (-zkr,). (4) 

Here Ek = m is the energy of a free particle, 

ek = {E^ + 2nQU EkY^"^ is the Bogolyubov energy of an 
excitation, and V is the normalization volume. The field 
operator of the force center can be represented in the 



form <l(r) — (l/Vl/) X]q ^xp (iqr), where Oq is the 
creation operator for the center. Then, using Eqs. (|^) 
and (^), Eq.(^) is transformed to 



k,q 



^—^ 6m 



where 
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qflq, (5) 



(6) 



The first term in the rhs of Eq. (^ , originating from the 
interparticle interaction, couples the motion of the force 
center with the excitation subsystem and is responsible 
for creating/annihilating excitations in the recombina- 
tion process. Considering this term as a small pertur- 
bation, we see that a single recombination event can be 
accompanied by the creation/annihilation of one excita- 
tion. Initially the momentum of the force center q = 0, 
and after the creation of the excitation with momentum 
k (annihilation of the excitation with momentum — k) 
the center acquires the momentum — k and the kinetic 
energy Ek/i. In a single recombination event, the prob- 
abilities of creating and annihilating the excitation with 
momentum k are given by 



w(iVk 



«;(^k ^ A^k - 1) 



1 |/tkP(l + ^k) 

V {ek + Ek/Zy 
1 |/^kP^k 

V [ek - Ek/if ' 



where A^k = [exp(efe/T) — are the equilibrium oc- 
cupation numbers for the excitations at a given temper- 
ature T . Then, for the rate constant of recombination 
accompanied by the creation of excitations we obtain 



-a 



<fik 



l + Nk 



{ek + Ek/^Y {ek-Ek/if 



(7) 



with a being the total (event) rate constant of recombi- 
nation. The first term in the rhs of Eq. (^ corresponds to 
spontaneous and stimulated creation of excitations, and 
the second term to their annihilation. 

For T <^ ^ the annihilation and stimulated emission 
of excitations can be omitted. One can put A^ = 0, and 
Eq.(0) yields 



(8) 



Even with a small value for the parameter (noa^) the 
large numerical factor in front of the expression (||) may 
imply that the creation of excitations in the course of 3- 
body recombination cannot be neglected. At the highest 
densities no « 3 x 10^^ cm~'^ of the MIT sodium experi- 
ment [|] Eq.(||) gives Qfox/a ~ 0.2. 

With increasing temperature, the role of annihilation 
of the excitations increases. However, our calculations 
from Eq. (|^) show that even at T ^ /i the annihilation and 
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stimulated emission of excitations give a small correction 
to Eq.(||). Only at T > lO/i the annihilation dominates 
over the emission, and aex (0) becomes negative. 

The most dramatic is the influence of created excita- 
tions on the loss rate of Bose-condensed atoms, which we 
will discuss for temperatures T < /i . For a high trap 
barrier single-particle excitations with energies ek ^ fJ- 
can be still trapped. These atoms act as "bullets" pene- 
trating the condensate. In a spherical trap this happens 
once per half of the oscillation period tt/u A charac- 
teristic time which a fast atom with velocity v^. spends in- 
side the condensate is ~ R/vk, where R = 
is the Thomas-Fermi radius of the condensate. Hence, 
the rate of elastic collisions of the fast atom with conden- 
sate atoms is ~ noavk{ujR/vk) ~ naucs, with cr = Sttc^ 
being the elastic cross section, and Cg = (/x/m)^/^ the 
sound velocity. In each elastic collision the fast atom 
transfers on average a half of its energy to the coUi- 
sional partner and removes it from the condensate. One 
has then to deal with two energetic atoms, and so on. 
This cascade process continues until the excitation en- 
ergy becomes of order the chemical potential /i. Ac- 
cordingly, the number of lost condensate atoms will be 
~ efc//J., and the characteristic time of the cascade pro- 
cess, T « 2{noaCs)^^ log (cfe/^). At realistic densities the 
time T is much smaller than the characteristic recombi- 
nation time Tr ~ (anp)^^. 

The behavior of the excitations produced in the cas- 
cade process depends on the ratio T//i. At T ^ /i 
their damping time strongly increases at energies well 
below jjL (the decay rate is at least much smaller than 
/i(noa^)^''^ ~ noacs fl^) and is likely to exceed the re- 
combination time Tr- Therefore, these excitations mostly 
remain undamped and no longer influence the number of 
atoms in the (partially destroyed) condensate. 

Thus, one has a non-equilibrium "boiling" Bose- 
condensed sample: High-energy single-particle excita- 
tions, created in the recombination process, initiate a 
significant destruction of the condensate and the forma- 
tion of a non-equilibrium non-condensed cloud. The cor- 
responding loss rate of condensate atoms, v — j Ln^cPr, 
is determined by the rate constant of recombination- 
induced production of excitations with energies Sk 3> /i, 
magnified by approximately efc//i: 



L = 



(Pk 



ek + Ek/S 



f-k 

—7, 



(9) 



where the numerical coefficient 7 ^ 1. A precise value of 
7 depends on a detailed behavior of damping rates of the 
excitations and, hence, on the trapping geometry. 

The generated non-condensed cloud has energy ~ /i 
per particle and occupies the volume which is of order 
the Thomas-Fermi volume of the condensate. Similarly 
to the condensate, this cloud decays due to 3-body re- 
combination and, in this respect, the quantity v describes 
extra losses of (condensate) atoms from the sample. 

The integral in Eq.(P) is divergent at high energies. 



and one should put an upper bound = Eb, where Eb 
is the trap barrier. The inequality Eb S> ^ justifies that 
Eq.(|) indeed gives the loss rate due to the production of 
high-energy excitations {sk ^ /i). From Eq.(||) we obtain 



1/2 
7- 



(10) 



As /i = nQ„iU, where no,„ is the maximum condensate 
density, the rate constant L is independent of the number 
of Bose-condensed atoms. 

The direct loss rate of Bose-condensed atoms due to 3- 
body recombination is i^o =/ 3an^(Pr, as three atoms dis- 
appear immediately in each recombination event. Then, 
using Eq.(|l0|) and the Thomas- Fermi density profile 
no(r) = nom(l — r^/i?^), we express the total loss rate 
of Bose-condensed atoms, i^t = ^o+J Ln^d^r, through vq: 



vt = va 



216 , 3,1/2 / Eb 
llV27r V M 



1/2- 
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(11) 



The situation is the same at T ~ /i, if the cascade 
production of excitations with energies e ^ fi makes the 
quasiparticle distribution strongly non-equilibrium and 
prevents the damping of these excitations caused by their 
interaction with each other and with the thermal cloud. 
The number of excitations produced in one cascade pro- 
cess is^ Eb/ fJ-, and the number of thermal quasiparticles 
with e ~ /i is ^ [ji/TiuoY . Thus, under the condition 
Eb > lJ,Nf^ one also has a non-equilibrium "boiling" Bose- 
condensed sample, and the loss rate of condensate atoms 
will be determined by Eq.(pT|). 

Our calculations assume the s-wave scattering limit 
k\a\ <C 1 p^ , and hence the maximum trap barrier for 
which they are valid is Eb « ?i^/2ma^. In the case of 
®^Rb, the triplet scattering length is a = 5.8 nm, and 
Eb = 75 ^K. Assuming 7 « 1, this gives Lt ~ 3Lo and 
shows a qualitative significance of our mechanism: the 
loss rate of Bose-condensed atoms is essentially magni- 
fied by the creation of high-energy excitations and their 
destructive influence on the condensate. To be more 
quantitative, one has to consider the kinetics of excita- 
tions produced in the sample by the initially high-energy 
(trapped) atom. 

In ongoing BEG experiments a characteristic tempera- 
ture of a Bose-condensed sample is in the range from 100 
nK to 1 iiK and, hence, the above estimated magnifica- 
tion of the loss rate of the condensate atoms (factor 3 for 
Eb ~ 75/iK) practically corresponds to switching off the 
evaporative cooling. With evaporative cooling on, the ra- 
tio Eb //i for temperatures smaller than fi is in practice 
ranging from 2 to 5. Then, at typical densities uq ~ lO^'' 
cm^'^ Eq.(pT|) only gives a 10% increase of the total loss 
rate of Bose-condensed atoms compared to Lq. To some 
extent this explains the recent experiments [pT| , where a 
strong increase of the 3-body losses in the condensate has 
been observed after switching off the evaporative cooling. 
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For T ^ II one can also think of the situation, where 
the cascade production of excitations with energies of or- 
der Id does not significantly destroy the equilibrium dis- 
tribution of quasiparticles in the sample. This should be 
the case ii Eb ^ I^N^. Then the damping of these ex- 
citations comes into play, continuously decreasing their 
energy and partially refilling the condensate. This damp- 
ing originates from (inelastic) scattering of a thermal ex- 
citation on a given excitation, which transfers them to 
the condensate particle and the thermal excitation with 
higher energy |15| , p^ . A characteristic damping rate is 
of order £:(noa'^y'^ , and even for the lowest excitations 
(e Hlo) it can be larger than the rate of recombination. 

Consequently, one can conclude that the energy of ex- 
citations produced in the recombination process is ther- 
malized in the gas. The Bose-condensed sample will be 
in quasiequilibrium, with continuously increasing tem- 
perature. This provides extra losses of Bose-condensed 
atoms. Due to refilling the condensate in the course of 
damping of the excitations, these losses will be smaller 
than the extra losses described by Eq.(|l0|) in the case of 
a non-equilibrium "boiling" condensate. 

The rate of energy transfer from the excitations, pro- 
duced in the recombination process, to the thermal cloud 
determines the increase of the internal energy U of the 
gas. One can write it as ?7 = / WnQ(Pr, where the quan- 
tity W is obtained in the same way as Eq.(Rh: 
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(2^ 



ek + Ek/3 



(12) 



Relying on Eq. ( p^ ) and the known expressions for U and 
the number of Bose-condensed atoms A^o as functions of 
T and the total number of particles (see |l^]), we have cal- 
culated the extra losses of condensate atoms \dNQ/dT\T, 
related to the increase of temperature. At initial density 
no ~ 10^^ cm~^ they do not exceed 10%. 

In conclusion, we have found that inelastic coUisional 
processes in Bose-Einstein condensates can be accompa- 
nied by the creation of elementary excitations. It is worth 
mentioning that this phenomenon is not related to BEC 
as itself. It originates from the presence of the mean- 
field interparticle interaction and will also occur in a non- 
condensed ultra-cold gas, as soon as the parameter na^ 
is not extremely small. We have revealed the influence of 
the production of high-energy excitations in the course 
of 3-body recombination on the loss rate of atoms from a 
trapped condensate. This effect is especially pronounced 
for a high trap barrier Eb, and it would be valuable to 
perform a systematic experimental investigation of the 
loss rate of condensed atoms as a function of Eb- 
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